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1.0 INTRODUCTION 


This User Manual documents the FORTRAN program RANDOM2. RANDOM2 is 
based on fracture mechanics using a probabilistic fatigue crack growth model. It predicts 
the random lifetime of an engine component to reach a given crack size (see Section 2.0, 
Theoretical Background). 

Included in this Manual are details regarding the theoretical background of 
RANDOM2, input data instructions and a sample problem illustrating the use of 
RANDOM2. Appendix A gives information on the physical quantities, their symbols, 
FORTRAN names, and both SI and U.S. Customary units. Appendix B includes 
photocopies of the actual computer printout corresponding to the sample problem. 
Appendices C and D detail the IMSL, Ver. 10 l , subroutines and functions called by 
RANDOM2 and a S AS/GRAPH 2 program that can be used to plot both the probability 
density function (p.d.f.) and the cumulative distribution function (c.d.f.). 
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2.0 THEORETICAL BACKGROUND 


Fatigue crack growth data are usually presented as cycles, N, to reach a particular 
crack length, a. The initial crack size is ai. It is generally accepted that under constant 
amplitude alternating stress, fatigue crack growth can be related to stress intensity through a 
first order differential equation. 3 

da/dN = C(AK) m (1) 

where C is a material parameter, m is a material property (often a constant) and AK is the 
stress intensity range. Stress intensity range is given by 

AK = YAo/ia 

where Y is a constant dependent upon component and crack geometry and Ao is the 
constant amplitude alternating stress. Therefore, equation (1) can be written as 

da/dN = C(YAo/ia ) m 


or. 


da/dN = C Y m Ac m 7t m/2 a my2 (2) 

Equation (2) can be integrated, from the initial crack length, ai, to the final crack 
length, af, to yield N, the number of cycles. The result is 


Thus, equation (3) gives the "cycles to reach a given crack length." 

Metallurgical evidence indicates that casting pores play a significant role in the high- 
cycle fatigue life of cast nickel base-superalloys, especially at high temperatures. 4 The 
location and size of these fatigue crack-initiating pores vary greatly from one aerospace 
propulsion system component to another. This accounts for the large variability in fatigue 
life and leads to consideration of fatigue crack growth as a random phenomenon. 

Fatigue life directly relates to casting pore size, and pore size can be used to 
determine initial crack size, a;. Thus, utilizing principles of both probabilistic analysis and 
fatigue crack growth, a quantitative probabilistic constitutive relationship between fatigue 
life and fracture mechanics parameters can be developed. Using the "randomized equation" 
approach, the fatigue crack growth model, given by equation (3) has the following form: 

N = f(C,m,Ac,a i ,a f ,Y) (4) 


_-m/2+l _-m/2+l 

f i 


-m/2+1 


(3) 
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or, in general, 

N = f(X.), i = 1 ,6, (5) 

where the Xi are the six independent variables in equations (3) and (4). Equation (3) is 
"randomized" by assuming the first four variables in equation (4) to be random. Assuming 
a small crack in a relatively large component leads to assuming Y = 1.0, a deterministic 
value. A deterministic final crack size was chosen since experimental evidence indicated 
that it was relatively unimportant 3 

Probabilistic analysis, via simulation, yields the distribution of the dependent 
random variable, cycles, N. A probability density function (p.d.f.) of cycles is generated 
using the maximum penalized likelihood method. Maximum penalized likelihood generates 
the p.d.f. estimate using the method of maximu m likelihood together with a penalty 
function to smooth it 3 
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3.0 INPUT DATA 


Data input for RANDOM2 is user friendly and easy to manipulate (see, for 
example, die file entided NORMAL.INP, in Section 4.0). Hie first five lines of input have 
the same format, namely 2E12.4, and the last two lines differ. The last two lines of input 
have the formats I3,2XJ3,2X,2E12.4,2X,I3 and 13, respectively. A brief line by line 
description is given along with an example for each line (Note: the ruler is to aid the user in 
formatting and is not a part of the input). A table listing the physical quantities, their units 
and symbols is given in Appendix A. 

1. Random Number Generator Seed, ISEED, and Sample Size, NTOT 
EXAMPLE: 

123456789012345678901234567890 
1 40 


2. Material Property, RMM 
EXAMPLE: 

123456789012345678901234567890 
28.0E-01 1.4E-01 


3. Initial Crack Size (Pore Diameter), RAI 
EXAMPLE: 

123456789012345678901234567890 
300.0E-06 45.0E-06 


4. Material Property, RCC 
EXAMPLE: 

123456789012345678901234567890 
2.20E-11 0.22E-11 


5. Stress Range, DELSIG 
EXAMPLE: 

123456789012345678901234567890 
6.2E+02 6.2E+01 
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6. The DESPL 1 parameters are NODE, INIT, ALPHA, EPS, MAXTT and are entered 
in that order as follows: 

EXAMPLE: 

1234567890123456789012345678901234567890 
21 0 50.0E-01 10.0E-05 30 

7. The DESPL parameter, IOPT, is entered as follows: 

EXAMPLE: 

1234567890 

2 
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4.0 SAMPLE PROBLEM FOR RANDOM2 


The objective of tins program is to predict the random lifetime, to reach a given 
crack size for an engine component . The theory is based on fracture mechanics, using a 
probabilistic fatigue crack growth model (see Section 2.0, Theoretical Background). 
RANDOM2 input parameters are given in Table Al.l. Note that the first four parameters 
are random. Their means and standard deviations are input by the user. The last two 
parameters, Af and Y, are deterministic and are fixed internally by the program. They are 
equal to the values shown in Table ALL 


Table Al.l RANDOM2 sample problem input (SI units) 


FORTRAN 

Name 

Distribution 

Type 

Mean 

Standard Deviation 
(V alue) (% of Mean) 

RMM 

normal 

28.0E-01 

1.4E-01 (5%) 

AI 

lognormal 

300.0E-06 

45.0E-06 (15%) 

RCC 

lognormal 

2.20E-11 

0.22E-11 (10%) 

DELSIG 

lognormal 

6.2E+02 

6.2E+01 (10%) 

AF 

N/A 

2.0E-03 

N/A 

YY 

N/A 

1.0 

N/A 


The input is entered in the following format in a file entitled NORMAL.INP. 


1234567890123456789012345678901234567890 


28.0E-01 
300.0E-06 
2.20E-11 
6.2E+02 
21 0 


40 

1.4E-01 

45.0E-06 

0.22E-11 

6.2E+01 

50.0E-01 


10.0E-05 


30 
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Execution of RANDOM2 (source code entitled NR2.FOR) produces an output file 
entitled RANDM22 giving intermediate results (see Appendix B). Execution also produces 
the plotfiles OUT1 and OUT2 (see Appendix B). These files are used to plot the X and Y 
axes of the probability density function (p.d.f.) and the cumulative distribution function 
(c.d.f.), respectively, generated by RANDOM2. The plots are drawn from the plotfiles by 
the SAS/GRAPH graphing program (see Appendix C). These plots for the sample 
problem are shown in Figures Al.l and A1.2. 

This same sample problem has been reported in Boyce and Chamis. 6 There, 
however, it utilized U.S. Customary units and an older version of RANDOM2 
(IMSL Version 9.2 subroutines). 



Fig. Al.l p.d.f. of log of mechanical cycles for fatigue crack 
growth model, using maximum penalized likelihood. 
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6.0 APPENDIX A 


PHYSICAL QUANTITIES, SYMBOLS, AND UNITS 

The physical quantities, their symbols, and units for the fatigue crack growth 
model are given in the following table. 


Table A1.2 Physical quantities, symbols, and units 
for fatigue crack growth model for RANDOM2 


Physical 

Quantity 

Theory 

Symbol 

FORTRAN 

Name 

Units 

SI 

U.S. 

Material Property 

m 

RMM 

m/cycle/MPa m 

in/cycle/ksi in 

Initial Crack Size 

A| 

RAI 

m 

in 

Material Property 

C 

RCC 

m/cycle 

in/cycle 

Alternating Stress 

Ac 

DELSIG 

MPa 

ksi 

Final Crack Size 

Af 

AF 

m 

in 

Geometry Dependent 

Y 

YY 

(dimensionless) 


Constant 
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7.0 APPENDIX B 


SAMPLE PROBLEM: SOURCE, INPUT AND OUTPUT FILES 
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8.0 APPENDIX C 


1. RNSET - 

2. RNNOR 

3. RNLNL - 

4. DESPL - 

5. GCDF - 


IMSL SUBROUTINE CALLS FROM RANDOM2 


Initializes a random seed for use in the IMSL random number generators. 

Generates pseudorandom numbers from a standard normal distribution 
using an inverse CDF method. 

Generates pseudorandom numbers from a lognormal distribution. 

Performs nonparametric probability density function estimation by the 
penalized likelihood method. 

Evaluates a general continuous cumulative distribution function given 
ordinates of the density. 
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9.0 APPENDIX D 


SAMPLE SAS/GRAPH (VER. 5.16) PROGRAM FOR RANDOM2 


data a; 

INFELE 'OUT1.CPR' FIRSTOBS=2;input x y; 

GOPTIONS DEVICE=HP7470; 
proc gplot; 

axisl label=(h=l f=simplex 'LOG OF CYCLES’) 
valuc=(h=l f=simplex); 
axis2 value=(h=l f=simplex) Iabel=none; 
plot y*x / haxis=axisl vaxis=axis2; 

TITLE H=1 A=90 F=SIMPLEX ’PROBABILITY DENSITY FUNCTION'; 
symbol i=spline v=square; 
dataB; 

INFILE ’OUT2.CPR’ FIRSTOBS=2;input x y; 
proc gplot; 

axisl label=(h=l f=simplex LOG OF CYCLES') 
value=(h=l f=simplex); 
axis2 value=(h=l f=simplex) label=none; 
plot y*x / haxis=axis 1 vaxis=axis2; 

TITLE H=1 A=90 F=SIMPLEX ’CUMULATIVE DISTRIBUTION FUNCTION’; 
symbol i=spline v=squarc; 
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